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The gluonic contribution to the nuclear modification factor Raa is investigated for central Au 
+ Au collisions at yfs — 200AGeV employing a perturbative QCD-based parton cascade including 
radiative processes. A flat quenching pattern is found up to transverse momenta of 30 GeV, which is 
slightly smaller compared with results from the Gyulassy-Levai-Vitev formalism. We demonstrate 
that the present microscopic transport description provides a challenging means of investigating 
both jet quenching and a strong buildup of elliptic flow in terms of the same standard perturbative 
QCD interactions. 



PACS numbers: 12.38.Mh,24.10.Lx,24.85.+p,25.75.-q 

The phenomenon of jet quenching [l[ has been one of 
the two striking discoveries at the Relativistic Heavy Ion 
Collider (RHIC) [2J, y] . It is widely believed to provide 
tomography of the quark-gluon plasma (QGP) created 
in ultrarelativistic heavy ion collisions. Furthermore, the 
strong elliptic flow, quantified by the Fourier parame- 
ter v 2 j explores the early bulk properties of the medium. 
Since ideal hydrodynamics can fairly reproduce the ob- 
served ^-dependence on centrality 4], viscosity in the 
QGP is believed to be small 0, [g]. Hence, the QGP 
behaves like a perfect fluid, which has been the second 
major discovery at RHIC. So far, it has not been possible 
to relate both phenomena by a common understanding 
of the underlying microscopic processes. 

Early attempts to describe the collective behavior by 
means of standard binary, elastic perturbative QCD 
(pQCD) collisions within a parton cascade have failed 
unless (unphysically) large cross sections are employed 
@, B Hi- With such cross sections, as we will also show 
in the following, the attenuation of jets would be far too 
strong, basically extinguishing all high-momentum par- 
tons. There exist several related theoretical schemes of 
the energy loss on the partonic level based on radiative 
pQCD interactions that are aimed at describing the ob- 
served amount of jet quenching HI U, 12, ljl 14. 15. lij]. 
In particular, the energy loss is found to be dominantly 
radiative. It is thus believed that different microscopic 
physics become manifest in perturbative interactions of 
jets with bulk particles on the one hand, and in seem- 
ingly much stronger interactions between medium con- 
stituents that give rise to the liquid-like bulk properties 
on the other hand. 

It has been recently demonstrated within BAMPS 
(A Boltzmann approach to multi parton scatterings) 
17, [ijj], a sophisticated kinetic pQCD based parton 
cascade including radiative contributions, that inelastic 
Bremsstrahlung processes can drive the system to equi- 
libration within a short time less than 1 fm/c. Moreover, 
pQCD Bremsstrahlung processes generate a sizeable de- 



gree of collective flow [19j and a small ratio of shear vis- 
cosity to entropy density r//s. Using parton- hadron du- 
ality the simulated elliptic flow coefficient v-i has been 
demonstrated to be in fair agreement with data employ- 
ing a coupling constant a s = 0.3 -j- 0.6. Some details 
on freeze-out |20|, hadronization and possible hadronic 
phase contributions have still to be addressed. For the 
created gluon plasma ry/s is in the range 0.08-^0.15, sup- 
porting the idea that the medium behaves like a nearly 
ideal fluid [H[2l|. 

The present study is motivated by the challenging 
question whether the employed treatment of pQCD inter- 
actions, including radiative processes, can simultaneously 
account for the quenching of high- momentum partons. 

Microscopic transport calculations provide a realistic 
way of understanding heavy ion collisions including the 
full dynamics of the system and allowing for non-thermal 
initial conditions. We employ BAMPS to investigate the 
evolution of gluon matter produced in heavy ion collisions 
at RHIC. Treating gluons as semi-classical and massless 
Boltzmann particles, elastic gg — > gg interactions are in- 
cluded via the leading order pQCD differential cross sec- 
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is used for inelastic gluon multiplication and annihilation 
processes. q_i_ and k^ denote the perpendicular com- 
ponents of the momentum transfer and of the radiated 
gluon momentum in the center of mass (CM) frame of 
the colliding particles, respectively. Gluon interactions 
are screened by a dynamically computed Debye mass 
m 2 D = d G na s J -f§^r^N c f, where d G = 16 is the gluon 
degeneracy factor for N c = 3 and / = f(p,x,t) is the 
local distribution function. 
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In the treatment of Bremsstrahlung processes the co- 
herent Landau, Pomeranchuk and Migdal (LPM) effect 
[22T ]. leads to a suppression of the emission rate for 
high energy particles. Such interference effects cannot 
be included in the microscopic transport calculations. 
Therefore, the theta function in JT]) is introduced to en- 
sure that only independent processes are considered, i.e. 
that only processes within the so called Bethc-Heitler 
regime are taken into consideration. Physically this im- 
plies that the formation time r of the emitted gluon 
must not exceed the mean free path of the parent par- 
ticle. The mean free path A g is given in the local co- 
moving frame of a computational spatial cell, which is 
not identical with the frame where the gluon is emitted 
purely transversal and thus r = l/k±. A Lorentz factor 
(I + f3 tanhy cos 9) has to be employed. j3 



7 = 



cosh y 



denotes the boost velocity from the comoving frame to 
the CM frame of the colliding particles, y is the rapid- 
ity of the emitted gluon measured from the CM frame 
and 9 is the angle between (3 and the axis of the collid- 
ing particles in the CM frame. For thermal energies the 
boost f} becomes small and the theta function reduces to 
0(fc^A g — cosh y). 
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FIG. 1: (Color online) Differential energy loss of a gluon jet 
in a static and thermal medium of gluons with T = 400 MeV. 
The contributions of the different pQCD processes to the total 
dE/dx are shown. Additionally dE/dx scaled down by a fac- 
tor of 4 is shown for a gluon jet interacting only via isotropic 
binary scatterings with a fixed cross section a — 10 mb. 

The cutoff in (fTJ) affects the overall value of the interac- 
tion probability for gg «-> ggg as well as the momentum 
shape of the spectrum of the radiated gluons. Fig. Q] 
shows the differential energy loss dE /dx of a high energy 
gluon jet in a static and thermal medium of gluons with 
temperature T — 400 MeV as a function of the jet energy 
E. Throughout this work we employ a fixed coupling 
constant a s — 0.3. The elastic energy loss from gg — > gg 
interactions is found to be only slowly rising with en- 
ergy to a level of dE/dx » 1 GeV//m at E = 40GeV 
and as expected the contribution from ggg — ► gg pro- 
cesses is negligible with dE/dx « 0.3 GeV/fm. Radiative 
gg — > ggg processes strongly dominate the energy loss, 



for instance resulting in a total dE/dx ps 29GeV/fm at 
E = 40 GeV, with dE/dx almost linearly increasing with 
the jet energy E. In comparison to the Arnold-Moore- 
Yaffe formalism [1J| a hard gluon with E = 40 GeV loses 
75 % of its energy about a factor of 2 faster [23| • 

The individual cross sections moderatly increase with 
the jet energy, yielding a total (crtot) ~ 3.5 mb at E = 
40 GeV 24| for a medium as characterized above. We 
note that a dependence of dE/dx on the distance the par- 
ton has propagated cannot be observed in our approach, 
since coherent effects are explicitely not included. 

For arguments to come, Fig. [T] also shows dE/dx 
(scaled down by a factor of 4) for a gluon jet interacting 
purely elastic with the medium at a fixed cross section 
of a = 10 mb and an isotropic angular distribution. The 
differential energy loss at this temperature is about a fac- 
tor of 3 -j- 4 larger than dE/dx from pQCD interactions, 
e.g. dE/dx\ a=10mh ps 94GeV/fm at E = 40 GeV. 

The initial gluon distribution for Au+Au collisions 
at 200 AGeV is chosen to be an ensemble of mini-jets 
with a lower momentum cut-off po = 1.4 GeV, produced 
in nucleon-nucleon collisions following a Glauber-model 
with a Wood-Saxon density profile and Gliick-Reya-Vogt 
parton distribution functions [25| with a A"-factor of 2. 
Free streaming is applied to regions where the energy 
density has dropped below a critical energy density of 
c c = 1 GeV/fm . With this setup it has been demon- 
strated lH, 3, 2(| that the experimental findings for the 
rapidity distribution of transverse energy and the flow 
parameter for various centralities are nicely reproduced. 

Due to the steeply falling momentum spectrum of the 
initial jets one would need to simulate an infeasibly high 
number of events in order to obtain sufficient statistics at 
high-pT- We have therefore developed a suitable weight- 
ing and reconstruction scheme. The underlying idea is 
to simulate a huge number of initial spectra, which is 
rather straightforward from the computational point of 
view, and to select such events for further simulations 
that contain high-px partons. The results then need to 
be appropriately weighted. For this we characterize each 
initial state according to X = max(pT,i); the maximum 
Pt in a given rapidity range. From each bin j in I a 
number of Nj events is simulated, the results are aver- 
aged within these bins and finally combined with appro- 
priate weights Pj. This procedure has been thoroughly 
tested [13] and, so far, allows for the investigation of ob- 
servables up to px ~ 30 GeV. In the following we use a 
bin-size of 1 GeV and select up to Nj = 120 events per 
bin for full simulation. 

Jet quenching is generally specified in terms of the nu- 
clear modification factor 
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which is the ratio of the particle yield in a A+A collision 
at given px and y to the yield in a p+p collision scaled by 
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FIG. 2: (Color online) Upper: Gluonic Raa at midrapid- 
ity (y e [—0.5, 0.5]) as extracted from simulations for central 
Au+Au collisions at 200 AGeV. The shaded area indicates 
the statistical error. For direct comparison the result from 
Wicks et al. fl6l ] for the gluonic contribution to Raa and ex- 
perimental results from PHENIX [1| for tt° and STAR [U 
for charged hadrons are shown. Lower: Gluonic Raa (trian- 
gles) for a scenario where all particles with pr > 8 GeV can 
only interact elastically with fixed a = 10 mb. 



the appropriate number of binary collisions. We directly 
compute Raa by taking the ratio of the final pr spectra 
to the initial mini-jet spectra. In this Letter we concen- 
trate on central (b = Ofm) collisions; Raa for non-central 
collisions can be straightforwardly studied and will be 



presented in an upcoming work 24| . Fig. [2] shows the 
result for the gluonic contribution to Raa, exhibiting a 
clear suppression of high-px gluon jets at a roughly con- 



stant level of R AA 



luons 



0.053. This constitutes the major 
result of the present study. 

The suppression of gluon jets is approximately a factor 
3-4-4 stronger than seen in experimental pion data. An 
excessive quenching, however, was to be expected since 
at present the simulation does not include quarks, which 
are bound to lose considerably less energy due to their 
color factor and dominate the initially produced jets from 
Pt ~ 20 GeV. As also discussed in [20| , the inclusion of 
light quark degrees of freedom would yield better agree- 
ment between BAMPS results and data on particle spec- 
tra, mean-pr and the pT-dependence of elliptic flow. In- 
deed, comparing with state of the art results from Wicks 
et al. II al. obtained within the GLV formalism Il3j, 




for the gluonic contribution to Raa (seen as the line in 
Fig. [2]) , which in their approach together with the quark 
contribution reproduces the experimental data, one finds 
better agreement. The gluonic suppression in their cal- 
culations is found to be Raa ~ 0. 07 -f- 0.09. They assume 
a gluon density per unit rapidity of dN g /dy = 1100 and 
the strong expansion of the fireball has to be included 
via an effective description. In the present microscopic 
treatment the full dynamics are selfconsistently included, 
where the initial gluon density is given by dN g /dy ~ 700, 
which then evolves to a final dN g /dy w 800. The inclu- 
sion of light quarks will provide important means of fur- 
ther verification and will be addressed in a future work. 
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FIG. 3: (Color online) Generation of elliptic flow in BAMPS 
at midrapidity for a noncentral collision with an impact 
parameter of b = 8.6 fm, obtained from calculations with 
a s — 0.3 and two freeze-out energy densities, e c = 
0.6 and 1.0 GeV fm -3 [2(j. For comparison the result em- 
ploying solely binary collisions with constant a — 10 mb 
(s c — 1.0 GeV fm -3 ) is shown. The grey-shaded band rep- 
resents the experimental value [30l ]. 

To demonstrate the importance of inelastic processes 
for explaining both elliptic flow and jet quenching within 
a combined framework, in Fig. [3] we compare the el- 
liptic flow for a scenario of only binary and isotropic 
collisions with a = 10 mb to the full results from [20(. 
Though the cross section is large, it still does not fully 
succeed in explaining the strong elliptic flow. A still 
larger one has to be (phenomenologically) employed in 
order to generate enough flow On the other hand, we 
also investigate a scenario in which jet-like gluons with 
Pt > 8 GeV can only interact elastically with the same 
fixed and isotropic cross section of a — 10 mb. In this 
calculation, the medium particles interact with the usual 
pQCD based matrix elements, including inelastic colli- 
sions. As already expected from Fig. [TJ the results shown 
in the lower panel of Fig. [2] explicitly demonstrate that 
already with a = 10 mb the suppression is too strong. 
So, while approaches using fixed and large cross sections 
might be able to reproduce thcrmalization times, angular 
correlations j2|| or elliptic flow 0, [27} , they would fail at 
simultaneously reproducing the correct nuclear modifica- 
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tion factor. 

The use of a full microscopic transport treatment also 
offers the possibility to investigate observables, that are 
not directly accessible in experiment. Fig. 2] shows the 
probability distribution of production points in the trans- 
verse plane of gluon jets with a final px > 10 GeV in a 
central (6 = fm) Au+Au collision using the full pQCD 
BAMPS version. Events in this illustration are rotated 
such that the jets emerge with their final pr parallel to 
the a;-axis and directed in positive ir-dircction. The spa- 
tial extent of the initial production points is roughly con- 
trolled by the parameter Ra in the Wood-Saxon distri- 
bution n A (r) = n /(l + eS r - RA ^ d ), with R A = 6.37 fm 
in our case. A bias towards production points near the 
surface is visible as would be naively expected for the 
strong jet suppression seen in Fig. [5] The escaping jets 
have been produced on the average 1.7 fm away from the 
surface and have undergone (N23) ~ 0.6 radiative pro- 
cesses ((N22) ~ 2.0 and (N 32 ) ~ 0.1). Still, at such high 
transverse momenta a few jets can indeed probe the in- 
terior of the reaction volume, corroborating the notion 
that the inclusion of the full medium dynamics is crucial 
to the investigation of jet quenching. 




FIG. 4: (Color online) Probability distribution of production 
points of jets with final pr > 10 GeV. Events are rotated such 
that jets in the final state travel in the positive ^-direction. 

For the first time a consistent and fully pQCD based 
microscopic parton transport description has been ap- 
plied as a common setup to both elliptic flow and jet 
quenching observed in RHIC experiments. As established 
previously, the gluon matter simulated in the parton cas- 
cade BAMPS exhibits a sizeable pressure build-up and 
a small ratio of shear viscosity to entropy, 77/ ' s. In this 
work we found that the suppression of high-pT gluon jets 
is roughly constant at R A l A ° ns ~ 0.053 over a large pr 
range, with inelastic gg ggg processes being the main 
cause of energy loss. The observed level of jet quench- 
ing is in reasonable agreement with recent analytic re- 
sults, though the suppression of gluonic jets appears to be 
somewhat stronger. Future studies including light quarks 
and a fragmentation scheme into the simulation will help 



to further confirm our findings that bulk properties of 
the medium and the energy loss of high-py partons can 
be described within a common transport approach. The 
investigation of elliptic flow at large transverse momenta 
will also provide important means of cross-checking. 

The authors would like to thank M. Gyulassy for stim- 
ulating discussions throughout this work. The simula- 
tions were performed at the Center for Scientific Com- 
puting of the Goethe University Frankfurt. 
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